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We compute the phase diagram of salt-free polyelectrolyte solutions using a modified Debye-Huckel 
Approach. We introduce the chain connectivity via the Random Phase Approximation with two 
important modifications. We modify the electrostatic potential at short distances to include a 
bound on the electrostatic attractions at the distance of closest approach between charges. This 
modification is shown to act as a hard core in the phase diagram of electrolyte solutions. We also 
introduce a cut-off on the integration of the modes of wave length smaller than the size over which 
the chains are strongly perturbed by the electrostatic interactions. This cut-off is shown to be 
essential to predict physical phase diagram in long chain solutions. 



I. INTRODUCTION 

Debye-Hiickcl (DH) theory derived nearly 80 years ago 
successfully explains many thermodynamic properties of 
dilute symmetric electrolyte solutions. 1 Recent computer 
simulations of hard-sphere ionic fluids 2 show that the DH 
approach for finite size ion solutions provides a rather 
accurate estimate of the critical temperature, although 
it underestimates the critical density. A conceptually 
simple modification of the DH approach, based on ac- 
counting for the interactions of ionic pairs introduced by 
Bjerrum 3 with residual ionic fluid, however, gives better 
estimates for the critical parameters. 

Though a large scientific effort has been made towards 
understanding the thermodynamics of simple electrolyte 
solutions, flexible charged chains solutions are not un- 
derstood. The Random Phase Approximation (RPA), a 
DH approach which includes the structure function of 
the chains, is widely used to describe polyelectrolytes 
solutions. 5 ' 6 RPA, however, provides unphysical coexist- 
ing curves in salt free polyelectrolyte solutions. 7 In par- 
ticular, RPA predicts that the critical density decreases 
with increasing the chains degrees of polymerization N, 
and it goes to zero as N goes to infinite. Moreover, the 
critical temperature increases rapidly as N increases sug- 
gesting that strongly charged chains are not water solu- 
ble in salt-free solutions. Recent computer simulations of 
charged chains reveal a critical density rather insensitive 
to the value of TV, and a critical temperature that does in- 
creases as N increases but not as rapidly as predicted by 
the RPA approach. 8 A one component plasma approach 
to polyelectrolytes solutions gives both a critical density 
and a critical temperature insensitive to the the degree 
of polymerization N. 9 The one component plasma ap- 
proach ignores the contribution from the fluctuations of 
the monomers. That is, the one component plasma as- 
sumes that the chains provide a constant charge density 
support over which the counterions fluctuate. It is diffi- 
cult to demonstrate and understand why and when it is 
possible to neglect the chain density fluctuations. More- 
over, the value of the critical temperature is strongly un- 



derestimated in the one component plasma approach to 
polyelectrolytes salt-free solutions. 

Here we describe a modified polyelectrolyte RPA 
model that predicts realistic values of the critical tem- 
perature, and a critical density rather insensitive to N 
in agreement with the simulations. The RPA approach 
described here includes two modifications. The first mod- 
ification is related to the cut-off in RPA. In salt free poly- 
electrolyte solution the electrostatic contribution to the 
free energy of charged Gaussian chains does not reduces 
to the DH limiting law. 5 Instead, a more strong electro- 
static contribution results in salt-free solutions that gen- 
erates instabilities in the free energy as the monomer con- 
centration goes to zero at room temperatures in strongly 
charged chains. 7 Scaling arguments 10-12 and computer 
simulations, 13,14 however, show that charged chains are 
stretched in monovalent ionic solutions on length scales 
smaller and of the order of the inverse screening length. 
Therefore, it is not possible to include the monomer 
density fluctuations within RPA, which assumes that 
the electrostatic interactions do not perturbed the chain 
statistics (a linear response theory). Here we introduce 
a concentration dependent cut off proportional to the in- 
versive screening length in the RPA electrostatic contri- 
bution to include only the contribution from the unper- 
turbed Gaussian chains. We also discuss other possibili- 
ties for the cut-off that yield similar modifications in the 
phase diagram. The second and less important modifi- 
cation involves including the contribution from the hard 
core of the ions and the monomers. In simple electrolyte 
solutions the hard core of the ions is included in the elec- 
trostatic contribution to the free energy to avoid spurious 
results as the concentration of ions increases. In this pa- 
per we modify the interaction potential and show that 
our modified potential describes the polyelectrolyte solu- 
tion at the same level of approximation as the modified 
DH for hard-sphere ionic fluids. 

In Section II we describe the modified RPA model for 
salt-free flexible strongly charged chains. In Section III 
we discuss the electrostatic interaction potential used in 
the model. In Section IV we state the problem of RPA 
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describing Gaussian chains and the required modifica- 
tions. In Section V we outline and discuss the results, 
and give some conclusions. 



III. INTERACTIONS BETWEEN CHARGES 

For dilute solutions of electrolytes the potential of in- 
teraction between charges can be taken in the form 



II. MODEL AND FREE ENERGY 

Let us consider a mixture of negatively charged poly- 
electrolyte chains and a residual ionic fluid. We assume 
that every chain consists of N monomers each carrying 
a negative charged — q and we denote the concentration 
of monomers as p- . Every ion in the system carries the 
charge +q and their total concentration is p + . Due to 
electroncutrality the following condition holds 



P- = P+ 

We write free energy of the system in the form 
F = F rcf + F cl 



(1) 



Here F re f is the free energy of the system without any 
electrostatic interactions and in the limit of good solvent 
it could be written as 
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where T and V is the temperature (in units of Boltzmann 
constant ks) and the volume of the system respectively. 
We also assumed that monomers of the chains and dis- 
sociated ions have the same sizes a, so we can define the 
reduced density p* as follows 



P* = a 3 (P-+P+) 



(4) 



The second term on the right hand side of Eq. (2) is the 
electrostatic contribution to the free energy. Using RPA 
approach we write it in the form 5 
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where the structure function of the chain g(k) has the 
following definition 
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Indices i and j in Eq. (6) run over all monomers of the 
chain and the average < ... > is taken over all possible 
chain confirmations. U(k) in Eq. (5) is the Fourier trans- 
form of the interaction potential between charges and the 
parameter JC is used to cut the integration over the large 
values of k. The purpose of the present work is to dis- 
cuss the RPA approximation, and determine appropriate 
forms of U(k) and JC. 
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where r is the distance between charges and e is the di- 
electric constant of the media. With this potential Eq. 
(5) calculated for a simple case of N = 1 and JC = oo 
reduces to — k 3 /12tt. That is, the electrostatic contribu- 
tion reduces to the limiting case {na < 1) of the DH 
expression for the free energy of hard-sphere ionic fluid, 



pDH .,3 

TV 12tt v ; 



(8) 



Here the inverse screening length k is given by k — 
(2) Anq 2 l B (p + + p_), the Bjerrum lenght l B = e 2 /eT and 
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(ln(l + .T)-a; 



(9) 



Clearly, Eq. (7) for U(r) can not be used in the RPA 
approach if na <~ 1. As the value of na increases the hard- 
sphere interactions between the charges become more and 
more important and they should be taken into account. 
For this purpose we introduce the following form of U(r) 



U(r) = — (1 -exp(r/o)) 
re 



(10) 



which suppresses electrostatic interactions on the scale of 
the size of the monomer. The Fourier transform of the 
above potential can be easily calculated 



U(k) = 



k 2 (l + k 2 a 2 ) 



(11) 



and for the electrostatic free energy given by Eq. (5) we 
get 
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The integral in Eq. (12) can be done analytically for the 
case of N = 1 and JC = oo which results in the following 
expression for the free energy 
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where 
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(13) 
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Another way to account for the hard sphere nature of 
the charges is to use the potential (7) but cut integration 
in (5) at K, = 1/a. This results in the following form of 
the free energy 



where 



7v = -sr*(~) 
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(15) 



(16) 



In Figure 1 we compare the results for r DH (x), t°°(x), 
and t k (x), given by Eqs. (9), (14), and (16) respec- 
tively. The curves presented in the figure show similar 
behavior as parameter x — na increases. In Figure 2 
we also show the phase diagrams obtained on the plane 
(p*,T* = a/Is) using different form of the electrostatic 
free energy given by Eqs. (8), (13), and (15). The fig- 
ure clearly shows that the phase diagram obtained using 
Eq. (13) is in reasonable agreement with Debye-Huckel 
theory where as the one obtained using Eq. (15) differs 
significantly. This fact leads us to the conclusion that 
the hard-core interactions between the charges could be 
correctly taken into account through the potential U(r) 
given by Eq. (10), but not through the large k cut-off 
integration parameter K = 1/a. However, in the next 
section we show that a different form of the cut-off /C 
is very important for the proper discription of polyelec- 
trolyte solutions. 



IV. POLYELECTROLYTE SOLUTIONS 

We first examine the dilute case of infinitely long poly- 
electrolyte chains. In order to calculate the integral (5) 
analytically we assume that the macromolecules obey 
Gaussian statistics. For g(k) defined by Eq. (6) we use 
the following approximation 



g(k) = l + 



N 



1 + Nb 2 k 2 /12 



~ 1 
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where b 2 is the mean-squared distance between the neigh- 
boring monomers of the chain. For simplicity we assume 
here that b — a. We also use U(r) in the form (7) valid 
for the dilute case ko«1 and we rewrite Eq. (12) with 
K, = oo in the form 
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Evaluation of the above integral leads to 
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The result F e i ~ k 3 / 2 for infinitely long Gaussian chains 
was previously obtained in a number of papers (for ex- 
ample sec Refs. 5, 17). 

With the use of Eq. (19) we take the second derivative 
of the free energy (2) with respect to the density p — 
p + + p- and in the limit of p — > we find 



dp 2 \TV J ~ 2p 
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The negative sign of the second derivative of the free en- 
ergy shows that dilute solution of infinitely long polyelcc- 
trolyte chains are unstable in good solvent conditions. 7 
In Figure 3 we show the phase diagrams on the plane 
(p* : T* — a/Is) obtained for different chain lengths N. 
The electrostatic contribution to the free energy used to 
obtain the diagrams is 
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and it is evaluated numerically. The figure clearly demon- 
strates that as N increases the region of instability be- 
comes broader and broader even at very dilute concen- 
trations. The same results can be obtain if we introduce 
a constant cut-of value JC* — JCa = 2ir. 

We suggest here that the unphysical behavior de- 
scribed above arises from the inappropriate contribution 
to the RPA electrostatic free energy at large k. At large 
k values the chains are perturbed due to the electrostatic 
interactions. RPA, however, is only valid if the interac- 
tions do not modify the chain conformations. There- 
fore, the cut-off has to be chosen such that we only 
include the electrostatic contributions that can be ac- 
counted by RPA. It is straight forward to determine the 
cut-off given that the major feature of polyelectrolyte 
solutions is the existence of a concentration dependent 
correlation length £. 10 On the scales smaller then £ the 
chain is stretched and on larger scales it obeys Gaussian 
statistics. For strongly charged polyelectrolytes £ can be 
estimated as 11 ' 12 



( p *y/2 
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In Fourier space scales smaller than £ are given by ka > 
27r/£. Since the contributions to the free energy on these 
scales can not be accounted within the RPA approach, 
we cut integral (22) 



K* = 2tt( P *) 



1/2 



(24) 



where 
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V. RESULTS, DISCUSSION AND CONCLUSIONS 

In Figure 4 we show the modified RPA phase diagrams 
obtained by evaluating the integral (22) up to ka = JC* 
with JC* given in (24). The modified RPA phase diagram 
has a much lower value of the critical temperature for 
large N than the standard RPA phase diagram Figure 3. 
The critical temperature initially increases as N increases 
in the modified RPA phase diagram, and it plateaus at 
very high N values at about ls/b ~ 7.3. Since most 
polymers have values of ls/b < 4.2 at room tempera- 
ture, salt-free solutions of polyelectrolytes with monova- 
lent counterions are water soluble at room temperatures 
in agreement with the experiments. 

Another important feature of the modified RPA phase 
diagram in Figure 4 is the prediction that the critical 
concentration remains nearly constant as N increases in 
agreement with the simulations. This is only true if the 
contributions from fluctuations of k > JC are completely 
ignored. That is, if we add the RPA electrostatic con- 
tribution from the 2n/a > k > JC in the free energy 
using the structure function for rod like units, though 
the critical temperature only increases slightly, the criti- 
cal concentration goes to zero as N increases. This sug- 
gests that the charge fluctuations at lower length scales 
do not contribute to the free energy or that they are 
strongly suppressed. The suppression of charge fluctua- 
tions on shorter length scales that £ can be explained 
by the fact that the electrostatic interactions arc not 
screened at these length scales. Therefore, the counteri- 
ons and monomers are strongly correlated. Indeed, these 
correlations stretched the chains. 

We argue above that in polyelectrolyte solutions the 
cut-of is a natural inverse length over which RPA breaks 
down. It is important to point out, however, that a con- 
centration dependent integration cut-off is also obtained 
in the one component plasma, 18 in which counterions 
fluctuate over a non-fluctuating charged medium. The 
physical reason and the concentration dependence of the 
cut-off in the one component plasma explained in Ref. 
18 is different than in polyelectrolytes. In a one compo- 
nent plasma the integration in the electrostatic free en- 
ergy given by the RPA approach should be carried over 
a finite number of the wave-vectors k in the same way 
as in the Debye theory of the specific heat of solids. 19 
Namely, the total number of degrees of freedom in the 
system 37V = 3pV should be equal to the total number 
of physically different modes with wave-vectors k within 
the spherical shell of radius JC. The number of modes is 
twice the number of the wave- vectors since each k has a 
sine a cosine mode. Therefore, one obtains 

™ (25) 

o 

which leads to 

JC* = (9tt 2 p) 1/3 (26) 



The above argument assumes that the important 
length scale is the average distance between ions d ~ 
2tt/JC ~ p 1 ^ 3 , and that this length scales imposes a 
"periodicity" even though the non fluctuating support 
for the ions is a solid with no structure. 

The distance between charges d <~ 1/p 1 / 3 is an impor- 
tant length scale to determine the break down of DH in 
simple electrolyte solutions. 20 Therefore, in simple elec- 
trolytes d is an important length scale as the system gets 
concentrated when the electrostatic energy between two 
charges separated by d, e 2 /(Aired) > K B T or ls/d > 1. 
Interestingly, the integration to a concentration depen- 
dent cut-off JC ~ p 1 / 3 versus integrating to the inverse 
hard core size of the ions JC ~ 1/a does not affect signif- 
icantly the RPA phase diagrams in simple electrolytes. 
Instead, in salt-free polyelectrolyte solutions a concentra- 
tion dependent cut-off is essential to recover reasonable 
values of the critical temperature as N increases and a 
N independent critical concentration. 

Though we argue that in salt-free polyelectrolyte solu- 
tions the natural length scale that determines the cut-off 
is £, we obtain very similar phase diagrams if we instead 
use d as shown in Figure 5. Notice that the importance 
of a concentration dependent cut-off is at dilute solu- 
tions where £ > d. Therefore, there should be counterion 
fluctuations at length scales d < I < £ that we did not 
include in our RPA approach. If we add a RPA contribu- 
tion from these counterion fluctuations (only the counte- 
rions) we get negligible corrections to the phase diagram. 
Our arguments are in agreement with the results in Ref. 
21 where a self-consistent polyelectrolyte RPA where the 
screening length from polyions is wavevector-dependent 
due to chain connectivity shows that at high wavevec- 
tors (short length scales), only counterions contribute to 
screening. At lower wavevectors, both polyions and coun- 
terions contribute to screening. 

In summary, we show that in salt-free solutions a wave 
vector dependent cut-off is essential to recover physical 
phase diagrams. The cut-off is determined by the ap- 
plicability limit of the RPA. In the presence of salt the 
electrostatic contribution to the free energy reduces to 
the well known DH limiting law. 5 In that case reasonable 
phase diagram are obtain without a concentration depen- 
dent cut-off, 22 although in principle one should only in- 
clude fluctuations from the monomers on the length scale 
where the RPA approach is valid (i.e., larger than £). 
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NAM) at Northwestern University. 
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FIG. 2. Phase diagrams of simple electrolyte solutions ob- 
tained using different forms of electrostatic free energy given 
by Eq. (8) (solid line), Eq. (13) (dash line), and Eq. (15) 
(dot line). 
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3. Phase diagrams of salt free polyelectrolyte so- 



lutions obtained within classical RPA approach for differ- 
ent lengths of the polymers N. From buttom to top 
N= 10,20,50,100, 1000. 
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FIG. 1. The plots for r(x) given by Eq. (9) (solid line), 

Eq. (14) (dash line), and Eq. (16) (dot line). 
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FIG. 4. Phase diagrams of salt free polyelectrolyte solu- 
tions obtained within modified RPA approach with K.* given 
by Eq. (24) for different lengths of the polymers N. From 
buttom to top N = 50, 100, 1000, 10000. 
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FIG. 5. Phase diagrams of salt free polyelectrolyte solu- 
tions obtained within modified RPA approach with /C* given 
by Eq. (26) for different lengths of the polymers N. From 
buttom to top N = 50, 100, 1000, 10000. 
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